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§ ! Abstract 
o 

We investigate the possible role that baryon number violating Yukawa interac- 
tions may take in the inclusive decay B — > X s ^y. The constraints, derived using the 
^ . experimental results of the CLEO collaboration, turn out, in many cases, to be more 



43 



5-H 



stringent than the existing bounds. 



> 

in 

The Standard Model (SM) of the strong, weak and electromagnetic interactions 
qq ■ is in very good agreement with almost all present experimental data, even though 

a few important predictions have not yet been tested. Still, most physicists would 
readily admit that the SM cannot be the final theory, both on aesthetic grounds as 
well as on account of certain well-founded technical objections. As a result, numerous 
Q_i! attempts have been and are being made in the quest of a more fundamental theory. 

Experimentally, there have been two main strategies to probe new physics. On the 
one hand, we attempt to directly produce, and observe, new particles at high energy 
colliders. On the other, we look for virtual effects of such particles and/or interactions 
in various low and intermediate energy processes. The decay b — > 57 is an excellent 
candidate for the latter option IH-fH . Experimentally, the branching ratio for the 



inclusive decay B — > X s 7 have been measured by CLEO [T5| and ALEPH [IB to be 

(1) 



f (3.15 ±0.93) x 10~ 4 (CLEO) 
BR(B X sl ) = \ 

{ (3.11 ± 1.52) x 1(T 4 (ALEPH) . 



The above are in good agreement with each other and with the SM prediction |I7| of 
BR(B -> A,7) = (3.29 ± 0.33) x 10" 4 . While a small window for the contribution of 
new physics does remain, this agreement can obviously be used to constrain deviations 
from the SM. 
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In this paper, we investigate the influence that a scalar diquark may have on the 
above decayQ. Diquarks abound in many grand unified theories (with or without 



super symmetry) and even in composite models [18]. While vector diquarks are con- 
strained to be superheavy^ (with masses of the scale of breaking of the additional 
gauge symmetry), no such restrictions apply to the masses of scalar diquarks. Con- 
sequently, such particles can be as light as the electroweak scale. For example, a 
diquark like behaviour can be found even in a low energy theory like the Minimal 
Super symmetric Standard Model (MSSM), albeit in the version with broken i?-parity. 



Diquark Type 


Coupling 


SU(3) C x SU(2) L x U(1) Y 


$1 




(6, 3, -|) 






(3, 3, -|) 


$ 3 


[hlfiQuYQij + h^(u Ri ) c d Rj ] $ 3 


(6, 1, -1) 


$4 


[hff(QLi) c Q Lj + ~hf{u Ri ) c d Rj ] $ 4 


(3, 1, -§) 


$5 


h\f (u Ri ) c u Rj $ 5 


(6, 1, -|) 


$ 6 


h[j (u Ri ) c u Rj $ 6 


(3, 1, -|) 


$ 7 


h$ (d Ri ) c d Rj $ 7 


(6, 1, |) 


$ 8 


h^(d Ri ) c d Rj ^ 8 


(3, 1, I) 



Table 1: Gauge Quantum Numbers and Yukawa Couplings of Scalar Diquarks (Q em = 

^3 + 1). 

A generic diquark is a scalar or vector particle that couples to a quark current 
with a net baryon number B = ±2/3. Clearly, this may transform as either a SU(3) C 
triplet or sextet. Concentrating on the scalars (for reasons mentioned above), the 
generic Yukawa term in the Lagrangian can be expressed as 

£y } = hfqXP LM ^ A + h.c, (2) 

where i, j denote quark flavours, A denotes the diquark type and Pl,r reflect the 
quark chirality. Standard Model gauge invariance demands that a scalar diquark 

1 A brief discussion on the sensitivity of the branching ratio B — > Xsj to scalar diquark-top 
contribution has been presented in jl7| . 

2 We do not consider the case of non-gauged vector diquarks as such theories are 
nonrcnormalizable . 
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transform either as a triplet or as a singlet under SU(2)i and that it have a U(l) 
hypercharge |Y"| = |, |, |. The full list of quantum numbers is presented in Table |l|. 

It is clear that the couplings h\f , h\f and hff must be symmetric under the 

exchange of i and j while hfff, h^f , hfj and h[j must be antisymmetric. For the 

other two, viz. h^f and h\f , there is no particular symmetry property. Note that 
the quantum numbers of $2,4,6 allow them to couple to a leptoquark (i.e. a quark- 
lepton) current as well. This implies that these particular diquarks could also mediate 
lepton-number (L) violating processes. Clearly, such leptoquark couplings need to be 
suppressed severely so as to prevent rapid proton decay. 

We make a brief interlude here to discuss the MSSM Whereas B and L are 
(accidentally) preserved in the SM (at least in the perturbative context), it is not 
so within the MSSM. Supersymmetry and gauge invariance, together with the field 
content, allow terms in the superpotential that violate either B or L [2(J. Catas- 
trophic rates for proton decay can be avoided though by imposing a global Z2 sym- 
metry [F2TJ under which the quark and lepton superfields change by a sign, while the 
Higgs superfields remain invariant. Representible as R = (— i) 3 - B - L + 2S ' ) where S is 
the spin of a field, this "i?-parity" is positive for the SM fields and negative for all 
the supersymmetric partners. However, while this symmetry is useful in prevent- 
ing phenomenologically unacceptable terms, it has no theoretical foundation and is 
entirely ad hoc in nature. Hence, it is of interest to examine the consequences of 
violating this symmetry, not in the least because it plays a crucial role in the search 
for supersymmetry. In our study, we shall restrict ourselves to the case where only 
the .B-violating terms are non-zero. Such scenarios can be motivated from a class of 
supersymmetric GUTs as well ||22|| . The corresponding terms in the superpotential 
can be parametrized as 

W 9 = K jk U R D R D R (3) 

where U R and D l R denote the right-handed up-quark and down-quark superfields 
respectively. The couplings \'^ k are antisymmetric under the exchange of the last two 
indices. The corresponding Lagrangian can then be written in terms of the component 
fields as: 

= Kjk (< dC A + utd*d c k + u*d c 3 dl) + h.c. (4) 

Thus, a single term in the superpotential corresponds to three different diquark in- 
teractions, namely two of type h^f and one of type h\f . 

The best direct bound on diquark type couplings comes from the analysis of dijet 
events by the CDF collaboration p3"l - Considering the process — * $a — ► <2Wj, an 



exclusion curve in the (m$ 4 , h^) plane can be obtained from this data. Two points 
need to be noted though. At a pp collider like the Tevatron, the uu and dd fluxes 
tend to be small and hence the bounds are relatively weak. This is even more true for 
quarks of the second or third generation (which are relevant for the couplings that we 
are interested in). Secondly, such an analysis needs to make assumptions regarding 
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the branching fraction of $a into quark pairs, a point that is of particular importance 
in the context of -R-parity violating supersymmetric models. 

There also exist some constraints derived from low energy processes. Third gener- 
ation couplings, for example, can be constrained from the precision electroweak data 
at LEP or, to an extent, by demanding perturbative unitarity to a high scale ||25|| . 
Couplings involving the first two generations, on the other hand, are constrained^ by 
the non-observance of neutron-antineutron oscillations or from an analysis of rare 
nucleon and meson decays [26, 27] . While many of these individual bounds are weak, 
certain of their products are much more severely constrained by the data on neutral 
meson mixing and CP-violation in the .fT-sector f2~8 |. It is our aim, in this article, 
to derive analogous but stronger bounds. 





Figure 1: Feynman diagrams that determine the one loop b — > s'y decay amplitude. 

Within the SM, the quark level transition b — ► sj is mediated, at the lowest order, 
by electromagnetic penguin diagrams shown in Fig. |I](a-c?). While only the top-quark 
diagrams have been shown, for consistency's sake, other charge 2/3 quarks should also 
be included. However, these contributions are negligible on two counts: (i) the small 
mixing angles and (ii) the corresponding loop integrals being suppressed to a great 
extent due to the smallness of the light quark masses. The matrix element for this 
process is then governed by the dipole operator: 

- ^V*V tb (^)C^ M (m w )sa, v F^[m b (l + l5 )]b. (5) 

The QCD corrections to this process are calculated via an operator product expansion 

3 Although many of these analyses have been done for the case of i?-parity violating models, 
clearly similar bounds would also apply to nonsupersymmetric diquark couplings as well. 
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based on the effective Hamiltonian 



AC s 

-^-v t :v tb J2c^)o^) 

V2 i= i 



(6) 



which is then evolved from the electroweak scale down to /i = /i& through renormal- 
isation group (RG) equations. A large correction owes itself to the chromomagnetic 
operator b — ► sG (G being a gluon) 



"''V, •., • , 9s , r ,SM 
- V ts V tb\ 



V2 



32w 



)C* M (m w )s a a^G^[m b (l + y 5 )]bp 



(7) 



which arises from the diagrams of Fig. §(a-c)), The Wilson coefficients Cj M (mw) 
and C% M {mw) can be evaluated perturbatively |2^-^T[ at the W scale where the 
matching conditions are imposed. The explicit expressions are 



C$ M (m w ) = x 



7 — 5x 



\r- x(3x — 2) 



C§ M (m w ) 
where x = mf/m^. 



x 



24(x - l) 3 

2 + 5a; — x 2 
8(x-l) 3 



+ 



4(x- l) 4 
3x 



In x 



A(x 



lnx 



(8) 



G 



W 



s b 



W 



G 



(b) 





s b 



t s 



G 3 



(<?) 



Figure 2: Feynman diagrams that determine the one loop b — > sG decay amplitude. 

The leading order results for the Wilson coefficients at the 5-meson scale, is 
given by 



G 7 (^t 



16/23 



C 7 SM (m w ) - 8 -C^(m w )[l - rT 2/23 ] + ^[1 - V~ 19/23 ] 
3 blS 



(9) 



5 



with r] = a s (mw) / oc s (fib) , calculated using the leading // dependence of a s , and the 
present world average value of the strong coupling constant viz. a s (mz) = 0.118 ± 
0.005. To this order, then, 



r(b -> sj) 



aG\m\ 
32tt 4 



V tb V t :C 7 SM (v 



(10) 



where a is the fine structure constant. As the above decay rate suffers from large 
uncertainties due to m& and the CKM matrix elements, it is prudent to normalise it 
against the measured semileptonic decay rate of the b quark 



r(6 



cez/p 



where z = m 2 c /m 2 and 



g{z) = i 



1927T 3 



Iz + 8z z 



K(z)g(z)\V ch \ 



12z 2 \nz 



(11) 



is the phase space factor. The analytic expression for k(z), the one loop QCD cor- 
rection to the semileptonic decay, can be found in Ref. j3"2"||. The explicit dependence 



on ml is thus removed, while the ratio of the CKM elements in the scaled decay rate 



viz. 



vsy t 



tb 



cb 



0.976 ±0.010. 



(12) 



is much better known than the individual elements. 

An updated next to leading order (NLO) analysis [|32| of the B — > X s 7 branching 
ratio with QED corrections has been presented in Ref. WJ\. Incorporating both the 



NLO QCD and the resummed QED corrections, the Wilson coefficient C^^fib) in SM 
can be expanded as 



cf(» b ) = cm + 



a s (Hb 
Air 



g Q 



(em) 



(13) 



For brevity's sake, we do not give here the expressions for C^^b) and C\ em> (iJ,b) 
as these can be found in Ref. [52] and Ref. |T7| respectively. The inclusion of the 
NLO and QED corrections in the b —>■ 57 decay rate has significantly reduced the 
large uncertainty present in the previous LO calculation. From the quark level b — > S7 
decay rate, it is possible to infer the B meson inclusive branching ratio BR(B — ► X^) 
by including the nonperturbative l/m& and l/m c corrections. These bound state 



i{em) , 
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corrections also have been taken into account in Ref. 

Having delineated the formalism, it now remains to calculate the additional con- 
tributions due to the possible presence of non-zero diquark couplings. At the one-loop 
level, the only new contributions to b —* 57 and b — > sG arise from the diagrams of 
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Fig. fj](e-/i) and Fig. @(d-<7) respectively^. In the generic case, apart from modifica- 
tions in the coefficients of 07,8, two additional operators arise. Denoted as 



7 



32vr 2 
9s 

327T 2 



sa^F^[m b {\ - 75 )]6 



(14) 



they differ from their SM counterparts in their chirality structure^. 

To keep the analysis simple, we shall assume that only one diquark multiplet is 
light and that all the fields within a multiplet are degenerate^]. With this simplifying 
assumption, the new contributions, at the electroweak scale, are given by 



C 7 D (m w ) 
C 7 D (m w ) 
C°{m w ) 
C£{m w ) 



No 
A 

Nc 
A 



{QzF.iy) + Q t F 3 (y)} hi*. + — {QtF 4 (y) - Q*F 2 {y)}r b l* s 
{Q^F^y) + Q t F 3 (y)}r b r* s + ^ {Q t F 4 {y) - Q*F 2 (y)} l b r* s 



A 



A' 1 
A' 1 

-AV2G F V*V tb ml 



{CsFM + C t F 3 (y)}l b l* s + ^ {C t F 4 (y) ~ C*F 2 (y)}r b l* s 

m b 

{C*F x {y) + C t F 3 (y)}r b r* s + ^ {C t F 4 (y) - C*F 2 (y)} l b r* s 

m b 



(15) 



where 



Fi(y) 
F 2 (y) 
F 3 (y) 



1 



12(2/ - I) 4 
2(2/- l) 3 
12(3/ - I) 4 

2(1 -y) 3 



6y 2 In y — 2y 3 — 3y 2 + 6y — 1 



y +2y\ny 



2 + 3y — Qy 2 + y 3 + 6y In y 



(16) 



3 - % + y 2 + 2 In jy 



with y = m 2 /m%. The color factor N c is —1 and 2 for triplet and 6 scalar respectively. 
Qt and Q$ are the charges of top quark and the diquark respectively. The color factors 
C t and C$ — for the diagrams of FigH/ and Fig.^? — are given in Table |2|. 

Note that, once again, we consider only such contributions, as involve the top 
quark. As is easy to ascertain from eq. fll5"D , for other quarks in the loop, the cor- 
responding integrals are too small to be of any consequence. Thus, any coupling to 

4 Clearly, to this order, none of $5.6 can mediate either of these processes and hence we shall not 
consider such fields any further. 

5 In all of these four operators, contributions proportional to the strange quark mass have been 
neglected. 

6 Largc splittings within a multiplet is disfavoured by the electroweak precision data. 
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Diquark Type 


c$ 


c t 


k 


h 


n 


r s 


$1 


-5/2 


1/2 












$ 2 


1/2 


1/2 












$ 3 


-5/2 


1/2 






u 33 


u 32 


$4 


1/2 


1/2 


h$/V2 


h$/V2 


n 33 


< l 32 



Table 2: Values of the coefficients C$, C t , h, h, "fb, t s for different types of Scalar 
Diquarks. 

the diquarks $ 7 8 , for example, would not be constrained to an appreciable degree by 
radiative b decays. In Table 0, we also display the relevant chiral Yukawa couplings 
(to the b- and s-quarks) for different choices of the diquark. 

In estimating the effects of scalar diqaurk couplings, it is useful to consider the 
ratios ]T7] £ 7 8 with 

_ C?{m w ) 

6 = 1 + cf M (m w y (17) 

and similarly for £ 8 . For the new operators 7 and 8 , we define, 

- _ C?(m w ) 

6 ~ C™{m w ) (18) 

and in an analogous fashion. With these definitions, the B — > branching ratio 
can be written as, 

BR{B^X sl ) = B 22 (6) + (e 7 +^)B 77 (6) + (es+i 2 8 )B 88 (6)+^B 27 (6) 

+ &B 28 {8) + {te s + £ 7 i 8 )B 78 {8). (19) 

In a parton level analysis, the photon would be monochromatic, with E y = E™ ax = 
mb/2. However, once the gluon Bremsstrahlung contribution is included, the photon 
spectrum becomes nontrivial and, for experimental purposes, one needs to make an 
explicit demand on the photon energy, namely 

E 1 >{1- S)E™ ax , (20) 

where 8 is the fraction of the spectrum above the cut. The values of Bij(8) are listed in 
ref. 0] for different choices of the renormalisation scale /ib and the cut off parameter 



on the photon energy 8. As is well known, some ambiguities exist in the choice of 



S 



fib which should, typically, lie in the region m&/2 to 2mj. For our analysis, we used 
fit, = m& and 5 = 0.9. We have checked that other values of 5 do not change the 
bound significantly. 

Since the new physics becomes operative only above the electroweak scale, the 
additional contributions to the operators O7 and Og will only serve to change the 
Wilson coefficients at mw Of course, the additional operators 7 and 0$ would 
influence the RG equations for C7 and Cg as well. However, since we are primarily 
interested in small C^ 8 (mw), it is safe to neglect any term in the RG equations 
involving these coefficients. 




-1 -0.6 -0.2 0.2 0.6 1 -0.012 -0.008 -0.004 0.004 

Product Product 



Figure 3: The partial width for B — > Xg^y as a function of the product of the diquark 
and/or R parity violating couplings for a fixed diquark mass of 100 GeV, The shaded 
region represents the la limits of the experimentally observed value. The curves for 
/133 /132 and h^h^ are identical to those for /133 /132 and h^hy} respectively. 

In the absence of an L-violating coupling, these diquarks clearly do not influence 
the semileptonic decay modes of the 5-meson. Thus, we may continue to normalize 
the radiative 6-decay against b — > cev e in order to avoid the severe dependence on 
mj. In Fig. |^, we plot the branching ratio BRiB — > X57) in presence of a diquark 
(multiplet) of mass lOOGeV. We continue to work under the assumption that only one 
pair of couplings are non-zero. Furthermore, we assume the said couplings to be realQ. 
That the curves should be parabolic in the product of the two couplings in question 
is obvious. To appreciate the fact that many of these curves have their minima lying 
on the SM value, one needs to consider the chirality structure of the corresponding 
diquark couplings (see Table 0). For example, combinations involving either of 
or ^,32 , imply that there is no left-handed coupling to the strange quark. From 
eqs.(|l"5D, it is then easy to see that the SM Wilson coefficients CV,8 remain unaffected. 

7 The extension to complex couplings is straightforward. The imaginary parts, however, can be 
better constrained from an analysis of the CP violating decay modes. 
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Consequently, the new contribution adds incoherently with the SM amplitude. For 
the rest of the combinations, though, the interference term is non-negligible leading 
to a shift in the minimum. Hence, unlike those for the first set of combinations (those 
involving ), the branching ratios corresponding to these sets are in agreement 
with the experimental numbers for two non- contiguous ranges of the product. 

In each of eqs.flTBD, the second term, whenever allowed, is clearly the dominant 
piece. This enhancement by the factor of m t /mb comes into play only when the 
diquark couplings to the bottom and the strange quarks have pieces with opposite 
chirality. In other words, for a diquark of the type $3 (or $4), the simultaneous 
presence of both the allowed types of couplings is severely constrained by the data 
on B — > Xg'-f. 



Products 
of Couplings 


Bounds from B — > X s 7 


la 


2a 


1.(1)1.(1) 

"33 a 32 
or h in h ™ 


[-1.0,-0.85] 
[-9.3 x 10- 2 ,5.8 x 10- 2 ] 


[-1.1,-0.74] 
[-0.2,0.12] 


00 oz 

or h$h$ 


[-0.16,0.25] 
[2.6,3.0] 


[-0.33,0.53] 
[2.3,3.2] 


1,(3)1(3) 
' 4 33 n 32 


[-9.2 x 10" 4 ,9.2 x 10" 4 ] 


[-1.4 x 10- 3 ,1.4 x 10~ 3 ] 


1(3)1.(3) 
lb 33 lh 32 


[-2.2 x 10~ 4 ,3.5 x 10~ 4 ] 
[3.3 x 10~ 3 ,3.9 x 10" 3 ] 


[-4.5 x 10" 4 ,7.3 x 10~ 4 ] 
[2.9 x 10" 3 ,4.1 x 10~ 3 ] 


1(3)1(3) 
rl 33 ,l 32 


[-0.12,0.12] 


[-0.18,0.18] 


1.(4)7(4) 
"■33 n 32 


[-2.5 x 10~ 3 ,2.5 x 10- 3 ] 


[-3.7 x 10- 3 ,3.7 x 10~ 3 ] 


1(4). (4) 
' 4 33 a 32 


[-1.1 x 10~ 2 ,-9.3 x 10" 3 ] 
[-9.4 x 10~ 4 ,5.9 x 10" 4 ] 


[-1.2 x 10- 2 ,-8.3 x 10~ 3 ] 
[-2.0 x 10^ 3 , 1.2 x 10~ 3 ] 


mm 

a 33 a 32 


[-0.35,0.35] 


[-0.51,0.51] 



Table 3: Limits on the Scalar Diquark Couplings for mq >i = 100 GeV. 

In Table |3], we capture the essence of Fig. |3| in the form of actual limits that can be 
set on such products of couplings, for a diquark mass of 100 GeV. Understandably, 
the 2a bounds are weaker than the la ones. Similarly, the color-sextet couplings 
are more severely constrained than the color-triplet ones. It should be noted that the 
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structure of the interaction terms h$ are the same as those for h on 



(3) 



Consequently, the 



fry \ 

bounds are exactly the same. A similar story obtains for h\j and h^ J . As discussed 
above, for a few of the products there are two non- contiguous bands allowed. For 
the combinations h^h^ and h^h^ , though, the second window (both at the la 
and 2a levels) lies beyond the perturbative limit and, hence, are phenomenologically 
uninteresting. 

As discussed earlier, h[f is analogous to the trilinear i?-parity violating coupling 



(4) 



\"j k . Thus the constraints on h^h^ are equivalent to those on the product A^A^. 
For each of these couplings, the best individual bound comes from the precision mea- 
surements at the Z pole pl],|55fl, and amounts to A(y 2 , A^-g < 0.50 at the la level. 
We thus do not do very well as far as this particular combination is concerned. This 
can be attributed to both the chirality and the color structure of the operator, each 
of which is "unfavourable" as far the b — > s'j decay is concerned. For most of the 
other combinations though, we do significantly better than the product of individual 
bounds 



24 
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Figure 4: The partial width for B — > X s ^ as a function of diquark mass, (a) For each 
curve, the associated product of diquark and/or R parity violating couplings is held 
to be 1, while all other couplings are set to be vanishingly small. The shaded region 
represents the la limits of the observed value. The curves for h^h^ and h^h^ are 



identical to those for h S3 h 



32 



and h^h^ 



respectively, (b) As in (a) ; but the non-zero 



products of diquark and/or R parity voilating couplings are held to be 0.005. 



In our effort to compare with the results available in the literature, we have, until 
now, held the diquark mass to be 100 GeV and varied the strength of its coupling. In 
reality, though a diquark is more likely to be somewhat heavier. For the sake of com- 
pleteness, we next investigate the dependence on the diquark mass (Figs. [|), while 
holding the product fixed. As is expected, the extra contribution falls off with m$. 
The fall-off is somewhat slower than mj 2 (see the expressions for Fi(y) in eq.([TB)) 
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and the effects persist till ~ 3 TeV. The different rates of fall-off are governed by 
the dominant Fi(y) in each case. The case for the combination h^h^ looks some- 
what nontrivial. However, the shape is just a consequence of accidental cancellations 
between various terms of eq . (|19|) . As the exact nature of these cancellations depend 
crucially on the value of the diquark couplings, not much should be read into the 
shape in general or the minimum in particular. 

The two dependences (m$ and coupling strength) that we have studied can be 
combined to rule out parts of the phase space. In Figs.||, we exhibit this for two 
particular combinations. In each, the shaded regions of the parameter space are in 
agreement with the experimental results at the designated level. For h^jh^, the 
second allowed region is beyond the perturbative limit. 




100 300 500 700 900 100 300 500 700 900 

m $1 (GeV) m^GeV) 

Figure 5: The region of the parameter space allowed by the data when all other cou- 
plings are set to zero. The lightly shaded area agrees with the data at la level, whereas 
the encompassing darker region agrees at 2a. 

In summary, we have studied the effects of the scalar diquark and/or R prity 
violating coupling to the branching ratio B — > X s 7. Among the possible new con- 
tributions, the scalar diquark mediated diagram yield promising effects. The precise 
measurement of this branching ratio at the upcoming B factories in near future and 
the reduction of theoretical uncertainty will improve the limits on the product com- 
bination of different scalar diquark and/or R parity violating couplings, we obtained. 
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